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Abstract 

Numerical simulations were conducted to 
investigate the performance characteristics of bleed 
through normal slots and its effect on the turbulent 
boundary layer development under zero and strong adverse 
pressure gradient caused by incident oblique shock. The 
solution to the compressible Navier-Stokes and k-s 
equations was obtained in a domain that includes the 
regions inside the bleed slot and plenum in addition to the 
external flow. The computational results demonstrate the 
interactions between the plenum, and bleed flow and the 
effect of incident shock on the boundary layer development 
downstream. The computed results agree with the 
experimentally measured pitot and static pressure 
distribution inside the slot. The bleed mass flow without 
incident shock was underpredicted over the range of 
plenum pressures. The computations predicted the 
measured increase in bleed mass flow with incident shock. 

Introduction 

Bleed is used in supersonic inlets to control the 
effects of flow separation associated with shock boundary 
layer interactions on ramps, cowls and side walls. In order 
to predict the effect of bleed on the boundary layer 
development and internal shock structure, few investigators 
[1-3] included ramp chambers, throat plenum and exit 
louvers in their supersonic inlet flow simulations. 
However, because of the enormous increase in 
computational time and grid generation complexity, most 
investigators simulated the global effects of bleed by 
changing the boundary conditions [4-7] and/or the 
turbulence models [8,9] in the bleed regions. The imposed 
mass flux in the bleed area was based on experimental 
measurements in some investigations [7] and on empirical 
correlation of the bleed discharge coefficient in others [4- 
6], Even when the turbulence model was modified for 
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mass removal at the porous wall [8] the effect of bleed on 
boundary layer development was not adequately predicted. 

Reproducing the experimentally measured boundary layer 
profiles over the porous wall even required the 
implementation of blowing rather than suction in the 
turbulence model [8]. Pavnter et al. [9] used a rough wall 
algebraic turbulence model to simulate the increase in the 
growth rate of boundary layer induced by bleed, by 
changing the roughness parameter in the numerical 
solutions to the Navier-Stokes equations until the boundary 
layer profiles matched the experimental measurements [10, 
1 1], They determined that roughness was a strong function 
of the upstream boundary layer mass flux removed, in the 
bleed rates between 3% and 15 %. However Paynter et al. 
pointed out that the experimental set-up whose data they 
used to develop their model, may have produced non 
uniformity in the mean velocity distribution and 
nonequilibirum boundary layer in the bleed region. In 
general, bleed models are restricted by the particular flow 
conditions and any inadequacies in the experimental data 
used in their model development. According to Hamed 
and Shang [12], the conclusions regarding the effects of 
bleed hole size and bleed location relative to the shock 
were not consistent among several experimental studies. 

Hamed, et al. [13-18] and Remlinger, Shih and 
Chyll [19-22] followed an alternative approach to 
investigating bleed effects. They conducted numerical 
simulations in which the viscous flow field was resolved 
inside the individual bleed holes and slots. This approach 
helps in understanding the important local phenomena that 
control the flow in the bleed regions, and how they are 
affected by the bleed configuration, and the external flow 
conditions. These type of investigations can and have been 
used in simple flow configurations to study the relative 
effects of bleed/slot angle, size and location. They 
revealed the presence of a separation bubble near the 
hole/slot entrance and a "bleed shock" that initiates inside 
the hole/slot and attaches to the downstream comer under 
certain conditions. The size of the separation bubble which 
controls the bleed discharge coefficient is strongly affected 
by bleed hole/slot slant angle and less by the free stream 
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conditions and plenum pressure. The boundary layer 
development downstream was strongly influenced by the 
bleed angle and the bleed shock strength. 

The purpose of present numerical study is to 
investigate the effects of external flow and plenum 
conditions on bleed performance. In the investigated 
configuration bleed is applied to a flat plate turbulent 
boundary layer through a normal slot. The numerical 
solution to the viscous flow inside the slot and plenum and 
over the plate surface were obtained for supersonic flow 
with and without impinging shock at the slot's upstream 
comer. The flow computations were conducted over a 
range of plenum pressures up to choked bleed conditions. 
The results, which are compared with experimental data, 
show the flow characteristics inside the bleed slot and the 
effect of bleed on the boundary layer development 
downstream with and without incident shock. 

Flow Configuration and Computational Details 

In the investigated flow configuration shown 
schematically in Fig. 1, bleed is applied through a normal 
slot to the turbulent boundary layer on a flat plate. The 
bleed slot whose width is one centimeter and its depth 2.54 
cm is connected to a 41.4 cm x 63.5 cm plenum. Outside 
of the slot, the solution domain extended 12.2 cm above 
and 40.9 cm along the flat plate surface. The boundaries 
AB and CD were located at 28.2 cm upstream and 12.7 cm 
downstream of the slot's upstream comer to match the 
locations of experimental velocity profile measurements in 
reference [6], 

The flow computations were conducted using the 
PARC code [23] with the compressible k-s turbulence 
model of Chien [24] modified by Nichols [25] for 
compressibility effects. The change in bleed mass flow 
rate was achieved through changing the static pressure at 
the plenum's outflow boundary. Referring to Fig. 1, the 
upstream boundary conditions consisted of free stream 
Mach number, stagnation pressure and stagnation 
temperature and a computed boundary layer velocity 
profile for flat plate computations which was matched to 
the experimentally measured profile's displacement and 
momentum thickness. The pre-shock and post-shock 
conditions were specified along the upper boundary to 
locate the inviscid incident shock at the upstream comer of 
the slot and flat plate surface. Flow variables were 
extrapolated at the downstream boundary and at the 
plenum's outflow boundary where the static pressure was 
specified to control the bleed mass flow. The initial 
conditions were specified from the solution without bleed 
for a given free stream Mach number. Subsequently, the 
computed flow field with the bleed slot and plenum was 
used to initialize the solutions at higher bleed mass flow 
rates (lower plenum pressure settings). 


The computational grid shown in Fig. 2 includes a 
308x68 grid over the flat plate surface, 89x68 grid inside 
the bleed siot and 187x148 grid inside the plenum 
chamber. Variable grid spacing was used in both x and y 
directions for grid clustering around the bleed walls, 
plenum chamber walls and at the plate surface, with 

x . = 0.9446x10“^ cm and y . = 0.4755 x 10~ 3 cm 
mm J mm 

corresponding to y + =0.959 at the inflow turbulent 

boundary layer at 2.46 free stream Mach number. 

Results and Discussion 

Typical results of the computed flow inside the 
slot near choked conditions are presented in Fig. 3 and 
compared with the experimental results of Davis et al. [27] 
reproduced in Figures 4 for incident oblique shocks whose 
inviscid impact point coincides with the slot's upstream 
comer at x=0. The figures present the pitot and static 
pressure contours for a shock generation angle a = 8°, at 
2.46 free stream Mach number. The computations predict 
very well both the shape and magnitude of the pressure 
contours inside the slot, the location of flow separation, 
and the size of the separation bubbles on the slot walls. 
The computed Mach number and velocity vectors, which 
are shown in Fig. 5, indicate flow angles of 33° from the 
plate surface, and high Mach numbers of 1.8 in a triangular 
region at the center of the slot opening. Disagreements can 
be observed between the experimental and computational 
results in this region. Since the pitot probe was always 
parallel to the slot walls in the experiment, one should 
expect discrepancies in these regions where the probe is 
oriented at 67° angle to the Mach 1.8 flow. The sharp 
gradients in the pressure across "bleed shock” were also not 
resolved by the measurements. The experiment indicates 
reattachment of the flow on the slot's downstream wall near 
the exit, which was not predicted by the computations. 

Effect of Incident Shock 

The effect of incident shock at the slot's upstream 
comer is demonstrated by comparing the computational 
resuits with and without the shock. Figures 6 and 7 present 
the computed Mach number and pressure contours inside 
the slot and in the neighboring regions over the plate 
surface and inside the plenum. The Mach number contours 
inside the slot indicate that when there is no incident shock, 
the separated flow region along the slots upstream wall is 
much larger while the separated flow region over the slot’s 
downstream wall near the exit is much smaller. In 
addition, the Mach number contours indicate less flow 
turning into the slot when there is no incident shock. This 
coincides with a weaker expansion fan at the slot’s 
upstream comer as can be seen from the pressure contours. 
Also without incident shock, the "bleed shock” which 
originates inside the slot is quickly weakened outside 
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through interactions with a second expansion fan at- the 
slot's downstream comer. 

The maximum computed bleed mass flow with 
and without incident shock was 5.1% and 2.35% of the 
incoming boundary layer for the slot geometry (D/5 = 
0.38). The pressure and Mach number contours indicate a 
separation bubble on the plate surface upstream of the slot 
in the case of incident shock. This is confirmed by the skin 
friction coefficient distribution of Fig. 8. One can see from 
Fig. 8 that there is a large difference in the predicted 
friction coefficients with and without incident shock within 
two slot widths downstream of bleed. In order to 
understand the cause of this difference, the mass flux, flow 
angle, and Mach number distribution across the slot 
opening are presented in Figs. 9-11. In the case of no 
incident shock, the mass flows out of the slot in the 
subsonic region behind the "bleed shock" and the 
downstream slot wall. This flow turns around the slot 
comer and forms the boundary layer over the plate surface. 
In the case of incident shock the subsonic flow goes into 
the slot opening in the region behind the "bleed shock", 
and the flow over the plate surface downstream originates 
outside the slot. 

According to the velocity profiles of Fig. 12, the 
effect of incident shock is seen to extend over two to three 
slot widths upstream and downstream of the slot, where the 
distorted velocity profiles are less full near the plate. The 
distortion upstream is associated with the shock induced 
pressure gradient and flow separation. The shape of the 
profiles downstream is affected both by the local flow 
conditions above the slot and by the strength of the "bleed 
shock" at the slot opening. Since the incident shock 
reduces the flow velocity, the profiles downstream of the 
slot are less full near the wall. The "bleed shock", has 
smaller effect on the downstream velocity profiles in the 
case with no incident shock because it is quickly weakened 
by the expansion fan at the slot's downstream comer. 

Plenum Interactions 

A considerable part of the computational effort 
was consumed in modeling the flow field inside the 
plenum. The increased computational effort was not only 
associated with the number of grid points inside the 
plenum (50.6% of the total) but also with the slower 
convergence of the flow at low subsonic plenum speeds. 
In order to determine the nature of the plenum interactions 
with the bleed and external flow fields, the bleed mass flux 
was monitored at the slot opening (bleed inflow) and exit 
(bleed outflow) during the computations. Fluctuations 
persisted in the bleed mass at both inflow and outflow 
boundaries when the computations were conducted using 
local time stepping. The amplitude of the bleed mass flow 
fluctuations was higher at higher plenum pressures and 


diminished near choking. Subsequently one time accurate 
flow simulation was performed with a=8° incident shock 
at choked bleed conditions. The computational results 
indicated no fluctuations in the bleed mass flow at the slot's 
inflow boundary' in this case and much smaller fluctuations 
in the bleed mass flow at the outflow boundary compared 
to the local time stepping results. Additional numerical 
computations with local time stepping were then performed 
without the plenum. The bleed mass flow in these 
computations did not fluctuate and coincided with the time 
accurate bleed mass flow at the slot's inflow boundary. 

The discharge coefficient, Q, and the flow 
coefficient, F, are used to present the bleed mass flow data 
in nondimensional form. They are defined as follows: 

I 7 r»i 

Q = m h / AM-Z-y (-^_) 2 fr-» 

* * 4 RT, ' y+r 

‘inf * 

When the stagnation free stream value P is 

t iof 

used for the reference pressure, P„ the discharge 
coefficient, Q, represents the ratio between the bleed mass 
flow and the ideal mass flow that could pass into the slot at 
sonic conditions if the flow was to expand isentropically to 
fill the bleed area A b . On the other hand, the flow 
coefficient represents the ratio between the bleed mass 
flow and the mass flow that would pass through the bleed 
area at the flow conditions p e U e outside the boundary 
layer 

F = m b / A b p u U u 

The predicted variation in the discharge 
coefficient with plenum pressure is presented in Fig. 13 
with the experimental results of references [26] and [27], 
The figure shows the computed bleed mass flow without 
plenum, as well as the range of bleed mass flow 
fluctuations in the computations with plenum using local 
time stepping. The computations without plenum agree 
with the experimental data in the case of incident shock, 
but under-predict the bleed mass flow without shock over 
the range of plenum pressures. The much higher discharge 
coefficient in the case of incident oblique shock can be 
caused by the increased static pressure above the slot 
and/or the velocity distortion upstream. In order to 
separate the two effects, a bleed discharge coefficient 
based on the local static pressure, P^, is presented in Fig. 
14. The value of the local static pressure was taken as P^f 
for no incident shock and Pj, the inviscid static pressure 
behind the incident shock (Pj/Pjjjf = 1.7). This figure 
indicates that the difference in the computed discharge 
coefficient with and without incident shock becomes 
smaller based on this normalization. The remaining 
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difference can be attributed to the incoming velocity 
profile distortion by the incident shock. 

Conclusions 

Numerical simulations were conducted for 
supersonic flow at Mach 2.46, with bleed applied through a 
normal slot with and without incident shock at the slot's 
upstream comer. The computed results agree with the 
experimental measurements inside the slot in the case of 
incident shock. The computations indicate higher flow 
angles at the slot opening and a narrower separation bubble 
on the slot’s upstream wall in the case of incident shock 
which leads to higher discharge coefficient. The "bleed 
shock" which originates inside the slot and extends as -an 
oblique shock outside is quickly weakened in the case with 
no incident shock by the interactions with an expansion fan 
at the slot's downstream comer. 

The plenum caused fluctuations in the bleed mass 
flux at the slot exit but time accurate computations at 
choked conditions indicated that these disturbances did not 
affect the flow and mass flux at the slot opening. The 
fluctuations were associated with the lateral movement of 
the bleed flow between the two separation bubbles on the 
slot walls. The computed bleed mass flow without plenum 
coincided with the predicted steady value at the slot inlet in 
the time accurate numerical simulations with plenum. The 
experimental data for bleed mass flow agreed with the 
computations in the case of incident shock but were higher 
with no incident shock. 
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Fig. 1 . Schematic of Computational Domain. 


Fig. 2. Computational Grid Near Bleed Slot. 
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Fig. 4. Pitot and Static Pressure Contours Inside Slot 
(Experimental Results [27]) M = 2.46, a = 8°. 
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(a) Mach Number Contours (b) Velocity Vectors 

Fig. 5. Velocity Field Inside Slot 



(a) With Incident Shock 


Fig. 6. Comparison of the Mach Number Contours With 
and Without Incident Shock 
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(b) Without Incident Shock 

Fig. 6. Comparison of the Mach Number Contours With 
and Without Incident Shock 



(a) With Incident Shock (b) Without Incident Shock 

Fig. 7. Comparison of the Pressure Contours With and 
Without Incident Shock 
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Fig. 12. Boundary Layer Velocity Profiles 
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